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We study a simple two-dimensional model for motion of an elastic fila- 
ment subject to internally generated stresses and show that wave-like propa- 
gating shapes which can propel the filament can be induced by a self-organized 
mechanism via a dynamic instability. The resulting patterns of motion do not 
depend on the microscopic mechanism of the instability but only of the fila- 
ment rigidity and hydrodynamic friction. Our results suggest that simplified 
systems, consisting only of molecular motors and filaments could be able to 
show beating motion and self-propulsion. 
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Cilia and flagella are hair-like appendages of many cells which generate motion and are 
used for self-propulsion and to stir the surrounding fluid. They all share the characteristic 
architecture of their core structure, the axoneme, a common structural motive that was 
developed early in evolution. It is characterized by nine parallel pairs of microtubules, 
which are long and rigid protein filaments, that are arranged in a circular fashion together 
with a large number of dynein molecular motors |3J . In the presence of a fuel which is ATP, 
the dynein motors attached to the microtubules generate relative forces while acting on 
neighboring microtubules. The resulting internal stresses induce the bending and a wave- 
like motion of the axoneme. 

These biological systems are complex, they consist of a large number of different compo- 
nents and various patterns of motion have been observed. Attempts to model their behavior 
are either based on the assumption that some unknown control system generates oscillatory 
motor activity || or that a self-organized mechanism is at work P,|J. Generically, the latter 
involves a dynamical instability. Theoretical studies of simple models for collective action 
of molecular motors have demonstrated the possibility of such instabilities [H,|fj7]]. Several 
examples of oscillatory motion of biological many-motor systems are known. Recently, it 
was suggested that spontaneous oscillations observed in muscles could be a property of the 
motor- filament system alone ^M- This idea is supported by the fact that the oscillations 
continue to exist after all regulatory systems are removed || but also by the observation 
that an in vitro motor-filament system shows the signature of a dynamic transition ||. 
Furthermore, the observations that flagellar dyneins are able to generate oscillatory motion 



on microtubules | 1Q| , and that isolated and de-membranated flagella in solution containing 



ATP above a threshold concentration swim with a simple wave-like motion |I]J support the 
idea that basic types of flagellar beating could result from a dynamic instability. 

In this article, we introduce a simple two-dimensional model which reveals many physical 
aspects of the motion of an elastic filament driven by internal forces, that should be relevant 
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FIG. 1. Snapshots of wave- like patterns generated by a motor-induced Hopf-bifur cation calcu- 
lated for different boundary conditions (solid lines): (A) Clamped head, position and slope are 
fixed. (B) Fixed head, position is fixed only. (C) Free head subject to a viscous load. The broken 
lines represent earlier configurations. The arrows indicate the direction of wave propagation. 



for flagellar beating. Our approach is inspired by studies of semiflexible filaments subject to 
external forces |12V[16|, however, in our case all motion is induced by internal stresses. Our 



model consists of two incompressible but elastic filaments of length L arranged at constant 
distance a L and rigidly attached only at one end which we call the head. A large number 
of molecular motors and passive elements holding the filament pair together are assumed to 
generate a coarse-grained force per unit length / which is an internal stress acting in opposite 
directions on the two filaments and induces the bending of the filament pair. The dynamic 
equations of this model define patterns of beating motion resulting from the internal forces 
which are assumed to oscillate. More interestingly, we show that characteristic wave-like 
patterns which propagate along the filament are generated most naturally by a dynamic 
instability of the motor-filament system, see Fig. [I] for examples. As we show below, the 
qualitative shapes of these patterns do not depend on the microscopic mechanism of force 
generation but only on the elastic properties of the filaments and on hydrodynamic friction. 
We demonstrate that these patterns lead to self-propulsion of the system and calculate the 
velocity of motion. 

In order to define our model and to derive the dynamic equations, we start from the 
enthalpy functional 
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-(C(s)Y + f(s)A(s)+A(s)(d s ^ 



ds (1) 



where r(s) is a parameterization of the shape of the filament pair by the arclength s, k 
is the bending rigidity, and C = n ■ d 2 s r is the local curvature with the filament normal 
n. The internal force density / couples to the relative local displacement of the filament 
pair A(s) = a ^C{s')ds' ||12|| . In order to impose the constraint of local incompressibility 
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{d s r) 2 = 1, we have introduced the Lagrange multiplier A(s). The equation of motion can 
be written as 

dt r = -^-nn + -ttj- (2) 

where tt and Tin are projectors on the filament tangent and normal, and we assume lo- 
cal anisotropic friction with tangent and normal coefficients £ii and respectively. The 
Lagrange multiplier is determined by the condition {d s r)d t d s r — p~3j] . 

In order to keep the description simple, we consider small deformations of a filament 
parallel to the x-axis, r(s) = (s + u(s),h(s)), which we describe by an expansion in the 
transverse and longitudinal displacements h and u. To quadratic order in d x h{x) we can 
write 



G 



2 (d 2 x h) 2 + af(x)(d x h(x)-d x h(0)) 



dx , (3) 



where we use the x-coordinate as parameter. We first discuss transverse motion which 
for small deformations is independent of longitudinal forces |T7|] and satisfies the equation 



C±dth = —K&^h + ad x f together with two boundary conditions at the head with x = 
and two conditions at the tail for x = L. We assume a free tail which implies d 2 h(L) = 
and nd x h(L) = af(L). At the head, we distinguish three different shown in Fig. 

(A) clamped head with h(0) = and d x h(0) = 0; (B) fixed head with h(0) = and 
Kd 2 h{fd) = —a J f(x)dx; and (C) a viscous load at x = with friction coefficient rj for 
which the condition on h(0) in (B) is replaced by r]d t h(0) = af{0) — K,d^h(0). 

We demonstrate the basic properties of this model, by first assuming that an oscillating 
force density with constant amplitude is generated by some unspecified mechanism: f m (t) = 
Re(/oe Ja;< ). The total force density / acting on the filament pair is the sum of the force f m , 
internal dissipative forces and in general the forces of elastic elements which locally connect 
the filaments. Introducing the complex Fourier amplitude h, where h(x,t) = R,e(h(x)e' lwt ), 
we can express the total force density as / = xv + / , where v ~ iua(d x h(x) — d x h(0)) 
is the complex amplitude of the local sliding velocity v = d t A. The coefficient x = + 
K/iuj) describes a viscoelastic response of the material between the filaments with dissipation 
coefficient A and elastic modulus K. The oscillating state is characterized by 

Kd^h — a 2 ituxdlh + t;±iujh = . (4) 

The homogeneous active force /q only enters via boundary conditions. Eq. (f|) and bound- 
ary conditions represent an inhomogeneous linear system which is solved by h = Ae kx l L 
leading to four complex values of k. The corresponding coefficients A are adjusted to satisfy 
the boundary conditions which leads to a solution with an amplitude proportional to the 
internal force h(x) ~ Jq. We can distinguish two different regimes: (i) hydrodynamic friction 
dominates \x\ 2 <C n^±/ua 4 ; (ii) internal viscoelasticity dominates \x\ 2 3> ^±/uja A . We can 
neglect in Eq. (f|), x m case (i) an d £j_ m case (ii). Fig. shows examples of the amplitude 
H and the gradient of the phase of h(x) = H(x)e~ l ^ x ^ for x — and different boundary 
conditions as dashed lines. The corresponding time dependent solutions 

h(x, t) = H{x) cos(u>t — (f>{x)) (5) 
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are propagating wave-like shapes qualitatively similar to those shown in Fig. |]. The sign of 
the local propagation velocity v p = u/d x </> of the phase allows us to determine the direction 
of apparent wave propagation. 

We have thus developed the framework to calculate and analyze wave-propagating so- 
lutions of our model and can now study motion generated by the properties of the motor- 
filament system via a Hopf bifurcation. We assume that the material between the two 
filaments which contains both molecular motors and passive elements has properties which 
can be characterized on a coarse-grained level by a nonlinear history-dependent response 
function. We will study the instability of a non- moving solution h(x) = towards wave-like 
patterns. For this case it is sufficient to consider only small amplitudes, \d x h\ C 1 as de- 
scribed above. Furthermore, in this regime the local sliding velocity v is small and we can 
ignore nonlinearities in v and restrict ourselves to the frequency dependent linear response 
/ = xv- Here, we have set the artificially introduced force fo = and characterize both 
passive and active material properties by the complex response function x(o>, Q) which can 
e.g. be calculated explicitly for a simple model || or measured experimentally flS ]. The 



out-of-equilibrium nature of the system is characterized by the parameter f2 which can for 
example be identified with the ATP concentration. Note, that for an active system x can 
have unusual behaviors which formally correspond to a negative friction (Re(x) < 0) or a 
negative elastic response (Im(x) > 0). 

In the case fo = 0, Eq. (|J) and boundary conditions becomes a homogeneous linear 
system which always has the solution h(x) = and which can now be reinterpreted as 
an eigenvalue problem for x- Spontaneous motion corresponds to nontrivial solutions to 
this problem. A discrete set of such solutions hi exists, each hi corresponds to a complex 
eigenvalue x — Xii^), i = 1, 2 . . . oo which we order according to < fc+iO^OI- 

Consider now a system initially at equilibrium with Q = 0. If O is increased, an instability 
occurs as soon as a critical value Q c is reached for which x(w c , ^ c ) — Xi(u c ) for a frequency u c . 
In the vicinity of this point, the system develops for Q > Q c motion with this frequency and a 
shape characterized by the nontrivial solution hi(x). This scenario applies to a supercritical 
bifurcation. Nonlinear terms of the response function and nonlinear corrections to the simple 
Monge representation can become important for larger Q, or they could change the nature 
of the bifurcation to subcritical. Typically, the instability occurs for the smallest value x = 
Xi(u;) since larger \x\ require larger values of Q which correspond to more system activity. 
Note, that the resulting pattern of motion is independent of the microscopic mechanism 
which leads to the instability. It is sufficient that the active material is capable to generate 
the response x — Xi EH- 

Fig. [2] displays examples of the amplitude and the gradient of the phase of h\{x) for 
boundary conditions (A) and (B), snapshots of the corresponding motion are shown in Fig. 
[I]. The boundary conditions play an essential role in selecting different types of motion. 
Observing the sign of d x (f> which determines the direction of the phase velocity we find that 
for clamped head (A) the wave propagates from the tail towards the head while for case 
(B) it propagates in the opposite direction. The amplitude H(x) also differs significantly 
between cases (A) and (B), see Fig. 0. The case (C) of a free head with viscous load r\ is 
similar to case (B) and therefore not shown in Fig. 0. For this example, the qualitative 
properties of motion induced by the dynamic instability are the same as those of the system 
driven by a homogeneous force f , see Fig. ^|. In fact, for the parameters chosen, |xx| 2 <C 
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FIG. 2. (a) Amplitude H(x) (in arbitrary units) of the wave-like motion characterized by Eq. 
d||) as a function of the position x along the filament axis for boundary conditions (A) as defined in 
Fig. [j] and £j_u;L 4 /k = 2500. (b) same plot for boundary conditions (B). (c) Gradient d x cj) of the 
phase along the filament axis for the same systems. The solid lines correspond to motion induced 
by a Hopf-bifurcation for the smallest response coefficient xii the broken lines to motion induced 
by a homogeneous internal force and x = 0- 
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K £,±/{ ujaA ) and the corresponding solution is not far from the solution for x — 0. The case of 
homogeneous force /o is simple and allows us to explain the effect of boundary conditions. A 
homogeneous internal force /o can be rewritten as boundary terms in the expression of the 
energy: G ~ af h(L) — af h(0) — aLf d x h(0) + J L dx (d x ]h) 2 K,/2. Its action is equivalent to 
two opposite transverse forces afo acting at both ends together with a torque aLf applied 
at the head. In the case of a clamped head this apparent force and torque are suppressed and 
the system is driven by a virtual force at the tail, propagating the wave towards the head 
|T6| . If the head is not clamped, the virtual oscillating torque at the head can propagate a 
wave in the opposite direction. 

Can these beating patterns propel the filament? Time-reversal symmetry has to be 
broken, h(x,—t) ^ h(x,t), for propulsion to be possible [pOR . According to Eq. (H), this 
requirement is fulfilled since d x <p ^ 0. Because of the symmetry h(x,t) = —h(x,t + 
there can be no net motion in the transverse direction. In order to estimate longitudinal 
motion, we have to study the displacement u(x). To second order in d x h, we can write 

u(x) ~ u(0) - - f X (d x h) 2 dx' , (6) 



2 Jo 

indicating that the dynamics of u(x) is governed by the motion h(x, t). Note that u(x)—u(0) 
is small, but the filament displacement w(0) can become large. The longitudinal component 
fi of the hydrodynamic force density — (£j_nn + £\\tt) ■ d t f acting locally on the filament is 
given by /; ~ — C,\\)d x hd t h — £\\d t u(x, t) in our approximation. The velocity of motion V 
is the time-average of dtu(0) and follows from the condition that the total longitudinal force 
vanishes. If an isotropic viscous load is attached to the head, this condition is J L fidx + 
r]d t u(0) = and we find V = Vq/(1 + T)/£\\L) where 

V = - ( ^ - l) £ f L H{x) 2 d x <P dx (7) 



Aw 

is the no-load velocity. If the head is not permitted to move, the filament generates a force 
F = Vo£\\L at the head. Note, that for isotropic friction both V and F vanish. For a 



semiflexible rod-like filament £±/£|| — 2 [21] and the direction of motion is opposite to the 

direction of phase propagation. 

The parameters chosen in Figs. |I] and ^| correspond e.g. to L ~ 40/xm, a ~ 20nm, 
~ 2 • 10 _3 Ns/m 2 , k ~ 4 • 10~ 22 Nm 2 , which is the elastic modulus of about 20 Microtubules 



P2| , and a frequency uj/(2tt) ~ 30s _1 . For this choice, we find a critical value Xi — 
(— 10 + 20i)Ns/m 2 . We choose an amplitude of h% with maximal value H/L ~ 0.1. In 
this case, the maximal local sliding velocity is v ~ 8/xm/s. In axoneme, dynein motors are 
spaced every 24nm along the microtubules. Assuming that only one microtubule pair is 
active, we estimate that for this choice a force per motor of 4pN corresponds to the critical 
value / ~ \xi\ v - This is a typical force created by molecular motors. Larger forces could be 
necessary to generate beating with larger amplitudes. Our result suggests that in this case 
several microtubule pairs could be active at the same time thus allowing for smaller forces 
per motor. Using the motion h\ obtained for boundary conditions (C) with a viscous load 
7] = 5 • 10~ 8 Ns/m, we find a spontaneous velocity of lateral motion V ~ 40/im/s, which 
is significantly larger than local sliding velocities and not far from experimentally observed 
values. 
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We have demonstrated that a pair of elastic filaments held together by an active, force- 
generating material, can induce wave-like patterns by a dynamic instability of the system. 
This study is motivated by biological flagella such as those of sperms which use such motion 
for self-propulsion. Our model suggests that the boundary conditions imposed at the ends 
select the type of beating pattern. This could be tested by micro-manipulation experiments 
which apply external forces and torques at the ends of beating flagella. We have restricted 
our study to a two-dimensional system, small deformations and the linear regime of the 
instability. Nonlinearities could play an important role in the regime of rapid self-propulsion 
with large amplitudes. Furthermore, the possibility of torsional motion in three-dimensional 
systems could allow for new types of behavior. 

The observation that wave-like motion can be generated in a self-organized way raises 
the idea that sophisticated control mechanisms may have evolved after the development of 
the basic axonemal structure in order to fine-tune the system and to create more complex 
types of motion. This concept suggests that artificially constructed systems consisting only 
of motors and filaments could already undergo beating motion and self-propulsion. 

We thank R. Everaers for stimulating collaborations, M. Bornens and H. Delacroix for 
introducing us to cilia and flagella, A. Parmeggiani and C. Wiggins for useful discussions 
and A. Maggs for a critical reading of the manuscript. 
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